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" \ ■ Abstract 

We suggest a class of generally covariant ghost-free nonlocal gravity models generating de Sitter 
c7^ , or Anti-de Sitter background with an arbitrary value of the effective cosmological constant and 

Cn ■ featuring a mechanism of dark matter simulation. These models interpolate between the general 

$H ' relativistic phase on a flat spacetime background and their strongly coupled infrared (Anti)-de 

Cm Sitter phase with two propagating massless graviton modes. 

<^ ■ 
^ . Introduction. It is widely recognized that one of the main challenges of modern physics is 

'^ ■ the problem of dark energy (DE) - the mechanism which is supposed to explain observable cosmic 

acceleration [T]. Numerous efforts aimed to reconcile evidences for this phenomenon with gravity 

theory ([51 [31 HJ [51 |B] , etc.) suffer from the fine tuning problem associated with the hierarchy of the 

cosmic acceleration scale vs the fundamental Planck scale. Though this problem served as a motivation 

to go beyond the simplest appropriate modification of general relativity (GR) - explicit cosmological 

j^ term, in this or that way it is creeping into almost all models of DE. Modulo certain exceptions 7 , 

most of them in fact look as a sophisticated way to incorporate into their action in addition to the 

Planck scale the horizon scale (whether it is a graviton mass of massive gravity |S] , multi-dimensional 

Planck mass in braneworld theories or the DGP scale in brane induced gravity models W, etc.). 

To circumvent this difficulty one could adopt another, perhaps more promising, line of reasoning. 

»y-\ ■ If a concrete fixed scale incorporated in the model is not satisfactory, then one could look for a model 

V^ . that admits cosmic acceleration scenario with an arbitrary scale. Then its concrete value compatible 

^T ■ with observations should arise dynamically by the analogue of symmetry breaking to be considered 

separately. Even this very unassuming approach is full of difficulties, because modified gravity models 

t~^ ■ featuring this property (like unimodular gravity 0, /(i?)-gravity [3], Lorentz breaking theory [3], 

^^ . etc.) generally violate some of its conventional symmetries and have additional degrees of freedom 

which might lead to ghost instabilities and make the theory inconsistent. This problem is central to 

numerous attempts to modify Einstein theory especially under the requirement to preserve its general 

covariance, and it will be a central question of this paper. 

►v> , Here we suggest a nonlocal infrared modification of the metric sector of the theory, which is likely 

^ ' to implement the above approach. It will be based on the realization of the old idea of a scale- 

5^ , dependent gravitational coupling - nonlocal Newton constant [TOl [HI [H] - and will amount to the 

construction of the class of diffeomorphism invariant, ghost- free models compatible with the GR limit 

and generating the de Sitter (dS) or anti-de Sitter (AdS) background with an arbitrary value of its 

effective cosmological constant A. In addition to fine-tuning argumentation of the above type, the 

driving force of our approach will be the understanding of the fact that, to resolve such issues of DE 

as cosmic coincidence problem, this scale cannot be encoded in the fundamental or effective action of 

the theory (like, for instance, explicit A-term, massive graviton or i? -I- R^ /A models [13]), but rather 

should arise dynamically by the analogue of symmetry breaking (see, for example, |7]). Moreover, as 

a bonus for the construction of the ghost-free cosmic acceleration we will also get in our model a new 

mechanism of dark matter (DM) simulation. 

Flat-space background setup. Here we begin our search for a nonlocal modification of the 
Einstein theory within the concept of the effective scale-dependent gravitational constant. At a qual- 
itative level this concept was introduced in [10] as an implementation of the idea that the effective 
cosmological term in modern cosmology is very small not because the vacuum energy of quantum fields 



(N 
> 



is so small, but rather because it gravitates too little. This degravitation is possible if the effective 
gravitational coupling constant depends on the momentum scale and becomes small for fields nearly 
homogeneous at the horizon scale. Naive replacement of the Newton constant by a nonlocal operator 
suggested in [1^ violates diffeomorphism invariance, but this procedure can be done covariantly due 
to the following observation [llj . The Einstein action in the vicinity of flat-space background can be 
rewritten in the form 

Se = ^jdxg^'^ (-^^"^G^. +0[i?^j) , (1) 

where G^u — R^d — ^g^j^^R is the Einstein tensor and l/D is the Green's function of the covariant 
d'Alembertian acting on a symmetric tensoi|j. This expression is nothing but a generally covariant 
version of the quadratic part of the Einstein action in metric perturbations h^i, on a flat-space back- 
ground. When rewritten in terms of the Ricchi tensor i?^^ ~ VV/i -|- 0[h?] this expression becomes 
nonlocal but preserves diffeomorphism invariance to all orders of its curvature expansion. 

Thus, the idea of a nonlocal scale dependent Planck mass [10] can be realized as the replacement 
of Mp by a nonlocal operator - a function M^(n) of D, 






MlRt^''- G^, ^ i?^"-^^ G^,, (2) 



which would realize this idea at least within the lowest order of the covariant curvature expansion. 
This modification put forward in |10illl| did not, however, find interesting applications because it has 
left unanswered a critical question - is this construction free of ghost instabilities for any nontrivial 
choice of M^(n)? Here we try to fill up this omission and put some constraints on A/^(n). 

To begin with, if we adopt the above strategy, then the search for M^(n) should be encompassed 
by the correspondence principle. According to it nonlocal terms of the action should form a correction 
to the Einstein Lagrangian arising via the replacement R ^ R + R^^'^ F {\3)G ^u ■ The nonlocal form 
factor of this correction F{p) should be small in the GR domain, but it must considerably modify 
dynamics at the DE scale. Motivated by customary spectral representations for nonlocal quantities like 
F{\I\) — J dm? a{m?) / {m? — U\) we might try the following ansatz, F{\3) — a/{m? — U\), corresponding 
to the situation when the spectral density a{m?) is sharply peaked around some m? . As we will see, 
for m? ^ this immediately leads to a serious difhculty. Schematically the inverse propagator of 
the theory - the kernel of the quadratic part of the action in metric perturbations h^^j - becomes 
^ — D -|- oL? /{m? — D) where the squared d'Alembertian D^ follows from four derivatives contained 
in the term bilinear in curvatures. Then its physical modes are given by the two roots of this expression 
- the solutions of the corresponding quadratic equation D = m?^ . In addition to the massless graviton 
with m?i = massive modes with rn^ = 0{m?) appear and contribute a set of ghosts which cannot 
be eradicated by gauge transformations (for the latter were already expended on the cancelation of 
ghosts in the massless sector - longitudinal and trace components of /i^i/). 

Therefore, only the case of m? = remains, and as a first step to the nonlocal gravity we will 
consider the action 

S=^ j dxg''^l^~R + aRt^'^^G,}i. (3) 

On the flat-space background this theory differs little from GR provided the dimensionless parameter a 
is small, |a| ^ 1. The upper bound on |a| should follow from post-Newtonian corrections in this model. 
The additional effect of a is a small renormalization of the effective Planck mass. Comparing the second 
term of ([3]) with ([T|) we have in the linearized theory the following relation S — — 2' I dx g^^^R + 



^We use sign conventions for the Einstein action in the Euclidean signature spacetime and curvature tensor conven- 



aO[/i^y], which ahows one to relate the constant M to Afp, 

AP = -^. (4) 

1 — a 

Treatment of nonlocality. At this point we have to discuss the treatment of nonlocality in ([3|) . 
In principle, handling fundamental theories with a nonlocal action is a sophisticated and very often 
an open issue, because their nonlocal equations of motion demand special care in setting boundary 
conditions. Contrary to local field theories subject to a clear Cauchy problem setup and local canonical 
commutation relations, nonlocal theories can have very ambiguous rules which are critical for physical 
predictions. In particular, the action ([3]) above requires specification of boundary conditions for 
the nonlocal Green's function l/D which will necessarily violate causality in variational equations of 
motion for this actiono 

To avoid these ambiguities and potential inconsistencies we will once and for all assume that our 
nonlocal action is not fundamental. Rather it is a certain approximation for the quantum effective 
action - the generating functional of one-particle irreducible diagrams - whose argument is the mean 
quantum field. This functional is necessarily nonlocal, and its nonlocality originates from quantum 
effects (by various mechanisms widely discussed in literature including fl4'). In this case boundary 
conditions for nonlocal operations are uniquely fixed by the choice of the initial (and/or final) quantum 
state, and manifest breakdown of causality in variational equations for this action is harmless under 
a proper treatment of their nonlocal terms. 

To begin with, this causality breakdown does not immediately signify inconsistency in the cal- 
culation of scattering amplitudes or in-out matrix elements. These amplitudes are determined by 
Feynman diagrammatic technique and do not have manifest retardation properties because they are 
not directly physically observable. Physically observable quantities like probabilities are bilinear com- 
binations of scattering amplitudes and can always be represented as expectation values ( in | O | in ) 
of certain quantum operators O in the initial quantum state |in). For example, the probability of 
transition from this state to some final state | fin ) , Pin->fin = ( in | fin ) ( fin | in ) — ( in | Pgn | in ) , is an 
expectation value of the projector Pgn = |fin)(fin| onto this final state. In contrast to in-out ma- 
trix elements these expectation values are subject to Schwinger-Keldysh diagrammatic technique [TS] 
which guarantees causality of ( in | 0{x) | in). This property can be formulated as a retarded response 
of this average to the variation of the classical external source J{y) coupled to the quantum fields in 
terms of which the observable 0{x) is built, 

5{\n\d{x)\\n) _^ . _. 



SJ{y) 



= 0, x"<y". (5) 



This property is also not manifest and turns out to be the consequence of locality and unitarity 
of the original fundamental field theory (achieved via a complex set of cancellations between non- 
local terms with chronological and anti-chronological boundary conditions). However, there exists a 
class of problems for which a retarded nature of effective equations of motion explicitly follows from 
their quantum effective action calculated in Euclidean spacetime |16| . This is a statement based on 
Schwinger-Keldysh technique TSj that for an appropriately defined initial quantum state |in) the ef- 
fective equations for the mean field g^^ — (in|g^^|in) originate from the Euclidean quantum effective 
action S = "SEucUdoaniff/^iy] by the following procedure |16)'^l . Calculate nonlocal S'Euciidoan[ff/iiy] and 
its variational derivative. In the Euclidean signature spacetime nonlocal quantities, relevant Green's 
functions and their variations are generally uniquely determined by their trivial (zero) boundary con- 
ditions at infinity, so that this variational derivative is unambiguous in Euclidean theory. Then make 



■^Taken literally with any choice of boundary conditions for l/D, the action Js} effectively symmetrizes the kernel of 
this Green's function, so that nonlocal terms in equations of motion never have retarded nature and, therefore, break 
causality. 

^We formulate this statement directly for the case of gravity theory with the expectation value of the metric field 
operator 3^1/(3;), though it is valid in a much wider context of a generic local field theory |16| . 



a transition to the Lorentzian signature and impose the retarded boundary conditions on the resulting 
nonlocal operators, 



r rr retarded 

"^Euclidean 



Sgt,i.{x) 



0. (6) 



These equations are causal {gfj,u{x) depending only on the field behavior in the past of the point x 
in full accordance with Eq.(|5])) and satisfy all local gauge and diffeomorphisni symmetries encoded in 
the original S'Euciidean[g^i.]Q 

We will assume that our model falls into the range of validity of this procedure, which implies a 
particular vacuum state |in) and the one- loop approximation (in which it was proven to the first order 
of perturbation theory in 17J and to all orders in [TB^). The extension of this range is likely to include 
multi-loop orders and the |in)-state on the (A)dS background considered below, for which this state 
apparently coincides with the Euclidean Bunch-Davies vacuum. 

Thus, the action Q is understood as the Euclidean one (this explains our sign choice in the 
Einstein term) with zero boundary conditions for l/D at infinity. It can be localized in terms of the 
auxiliary tensor field subject to the same Dirichlct boundary conditionf|j, and in the resulting local 
representation directly applied to the FRW cosmology. This shows that close to a certain moment Iq 
corresponding to the present epoch the model easily yields a (quasi) de Sitter point of the cosmological 
evolution [18], its Hubble factor H — a/ a and the equation of state parameter w = —\ — 2H/3H^ 
satisfying the relations w{to) = —I, w{to) = 0(1) x H{to) < 0, which make the model qualitatively 
compatible with the observable DE data. These preliminary estimates could have served as a starting 
point for a quantitative comparison with the DE scenario. However, a formal application of ([3]) to 
the FRW setup disregards nontrivial boundary conditions in cosmology. To see this, note that on the 
de Sitter background (which is a zeroth-order approximation for the cosmic acceleration scenario) the 
Ricci curvature i?^^ = Ay^^ is covariantly constant, and the nonlocal part of ([3|) is divergent, because 
g^y is a zero eigenvector of D. This means that the action ([3|) should be modified to circumvent this 
difficulty. 

Nonlocal gravity with a stable (A)dS background. We will regulate the action ^ by adding 
to the generally covariant D the matrix- valued potential term built of a generic combination of tensor 
structures linear in the curvature, 

S=^ [dx gi/2 (-R + a i?'^"^ ^1-. G^,) , (7) 
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Here we use the condensed notation for the Green's function of the operator D + P = D 5 A 
acting on a symmetric tensor field as 



a/3 



^ G = 



U + P 



1 



U + P 



G^a (9) 



and a, &, c, d and e represent arbitrary parameters to be restricted by the requirement of a stable 
(A)dS solution in the model. Of course, such a modification of the original action ([3]) leaves its linear 
approximation on a flat background intact, because it deals with 0[/i^^] -terms. 



* A similar treatment of a nonlocal action in |19| was very reservedly called the " integration by parts trick" needing 
justification from the Schwinger-Keldysh technique. However, this technique only provides the causality of effective 
equations, but does not guarantee the Euclidean- Lorentzian relation ||6ll. The latter is based, among other things, on 
the choice of the |in)-state. 

^This field formally carries ghosts, but this does not indicate physical instability because it never exists as a free field 
in the external lines of Feynman graphs. The actual particle content of the theory is determined in terms of the original 
metric field g^^ and indeed turns out to be ghost-free on the flat-space background, because the quadratic part of the 
action coincides with the Einstein's one. A similar mechanism excluding ghosts by boundary conditions was recently 
used in the conformal gravity model of [20) . 



Now the Green's function 1/(0 + P) acting on the Einstein and Ricci tensors in ([7]) is weh defined 
even for the (A)dS background with the covariantly constant i?^^ — h^g^u and Raf^pi, = ^{gapSf^u — 
Qaugp^i), for which 

{gcpg'"). (10) 
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(11) 

(12) 

so that Pg^,^ = F^^"'^ g^p = {A + 45) Ag^^. 

Let us show that under a certain restriction on parameters of P the model (O has (A)dS solution 
with an arbitrary value of the cosmological constant A. Indeed, introduce the local conformal variation 
with the parameter 5a = S(t{x), 6^ = J d'^x Scr gap S/Sgai3- It operates on various quantities in 
([7]) according to their conformal weights, S^g^^ = Sagfj^i,, S^R^j.^ — 0(V), S^R = —SaR + 0(V), 
Sa-P = —Sa P+0{V), Sa-O = 0(V), modulo the derivatives 0(V) acting on da'{x) and these quantities 
themselves. Therefore, the conformal variation of ([7]) on the (A)dS background reads 



5S_ 
1^ 



,«.s ^ ^""l-" + ° «'""^" "" ) - -'"'' (' + 7TS ) »"'■ C^' 



On this background the full metric variational derivative expresses via its trace part given by this 
variation, SS/Sg^^i, |(A)dS = j9^'^{^S/S(t) |(A)dS- Therefore, the equation of motion 5S/5g^i/ |(A)ds = 0; 
holds with an arbitrary value of A when 

a^-A-AB. (14) 

Note that the existence of the (A)dS solution with an arbitrary A is neither the result of the local 
Weyl invariance of the theory, nor even its global scale invariance. Rather this is a corollary of the 
relation ([T4)l which, in particular, guarantees the vanishing on-shell value of the action S'|(A)dS = 0- 
Thus, this solution is another vacuum - a direct analogue of the flat-space one. 

Another remarkable consequence of Eg. p^ is the stability of the (A)dS solution against ghost and 
tachyon excitations. In principle, the hope to eradicate ghosts and tachyons from the quadratic part 
of the action i5'(2) on the (A)dS background is based on the observation that in the DeWitt gauge, 

X^ = V^/i^'' - i V/i = 0, (15) 

S(2) contains only two contractions Kj^^ and h? (h = g'^'^h^i,), their nonlocal parts being given by 
h^"{\I\ + P)^^hf^i, and /i(n — aA)~^ft,0 As in the discussion above, the ghosts necessarily appear if 
these nonlocalities are nonvanishing in S(2) j because the dispersion equation for D becomes quadratic 
and generates doubled set of physical modes with D = rn^ {P playing the role of nonvanishing rra^ 
above). Therefore it is a priori possible to cancel these two nonlocalities and provide the right signs 
of the remaining local terms by an appropriate choice of five parameters in (|5]). 

Curious fact is that in the DeWitt gauge Si2\ depends only on the traceless part of ft.^^, /i^j/ = 
h^v — \ gtiiyh, and is very simple 

5(.)=-^/dW/^{/.-nV + (c-l)A/.^.+A^(c4)^M-^/.,.}, (16) 
^°«^8B(2i? + a)^^' ^^^^ 



^The form of the Green's function in the trace sector follows from the equation (D + P)g^vh = 9ti,v{^ — aK)h also 
based on II14I I. 



where Geff = l/87rM^ff is the effective gravitational constant vs the Newton constant Gjv = l/SiiMp. 
This is the main technical result of this paper (the detailed derivation of which will be reported in 
[TS]), and it leads to the single equation for the constant p^ . 

C^^-c-4e=| (18) 

and the positivity requirement for Goff. This requirement selects the range of the parameter _B, 
B < —a/2 and B > Q for a positive a, and even more interesting compact range of B for a negative 
a, 

0<B<--, a<0. (19) 



The action without gauge-fixing can be obtained from ([16)) by representing /i^i. in the DeWitt 
gauge as the projection of the non-gauged field, h^i,\^<^=o = /i^^ — 2 [V(^;^"j/(n -I- A)] Xa- The result 
reads 



5f2) = ^ /d^ 1/2 r 1 / 2 \ , l/^/'_g_2 \ 1 

'^' 167rGcff J -^14 V ay'^'^sV 3/ 2 



Au 



1 

4"^'^n + 2A 



^(i) n-L9A -^(i) r ' (^^) 



where, quite interestingly, the first three terms coincide with the quadratic part on the (A)dS back- 
ground of the Einstein- Hilbert action with the A-term and i?(i) — VpX^~ ^C^ + ^A)/! is the linearized 
Ricci scalar on the same background. 

The trace of the variational equation for ([20)) . according to ([6]), gives the homogeneous equation 
for i?(i) with a retarded nonlocality, which is equivalent to the local initial Cauchy problem □i?(i) = 
with zero initial data in the remote past [18]. Therefore i?(i)(a;) = 0, whence in the DeWitt gauge 
X'' = we have (D -I- 2A) h = Q. Then the residual gauge transformations /S.fh^y — 2V(^/i^) with 
the parameter /^ satisfying the equation (D + A)/^ = can be used to select two polarizations - 
non-ghost physical modes. In particular, the boundary conditions for h can be nullified, so that h 
identically vanishes and makes in view of the DeWitt gauge the propagating free modes transverse 
and traceless as in the Einstein theory with the A-term. 

In the presence of matter sources with a stress tensor T^j/ of a compact support the causal effective 
equations for retarded gravitational potentials become local in the DeWitt gauge, 

□ + ^AJ V + ^V^V,/i - -g^,,h = j^T^u- (21) 

Modulo the gauge transformation their solution takes the following form - the result of a careful 
commutation of covariant derivatives with (— □ + | A)^^, 

The tensor structure here differs from the GR analog T^^ — \g^^T , which for non-relativistic sources 
gives 0(1) correction. What is much more interesting, it yields an unexpected bonus in the form of the 
dark matter simulation - l/|a|-amplification of the gravitational attraction due to the replacement of 
the Newton gravitational constant Gat by ([TT)) . Goff ~ Gm l\ct\ with |a| <C 1. This is possible when 
\B\ ~ a for a positive a and necessarily happens in the case ()19p of a negative a, because the factor 
a/8B{2B + a) > l/|a| and 

1 -a 
Gcff > — i — j— Gat ^ Gat. (23) 



On the other hand, with \B\ ~ y/a in the case of a positive a both Newton and effective gravitational 
couphng constants can be of the same order of magnitude, Gn /GeS = 0(1) even for a <C 1, which 
together with ([^5]) leaves a large window for a possible strength of DM attraction relative to GR 
behavior. 

Conclusions. The theory ([7]) has two phases. For short distances corresponding to the range 
of wavelengths with VV ^ D ^ i? this is a GR phase on the zero curvature background with small 
0{a) X R/D corrections of higher orders in spacetime curvature (collectively denoted by R). This 
regime applies to galactic, Solar system and other small scale phenomena and is likely to pass all 
general relativistic tests for a sufficiently small a|3 

Another phase of the theory corresponds to the infrared wavelengths range VV ^ R in which a 
stable (A)dS background exists and the modified gravitational potential of matter sources is given 
by Eq.lgll). This equation is vahd for the perturbation range \SR^\ ~ |VV/i{^| < A and |/i(^| < 1 
equivalent to very small matter densities |T^| ^ MpA characteristic of galaxy, galaxy cluster and 
horizon scales for which DE and DM modification of gravity theory becomes important. Thus, nonlocal 
gravity interpolates between GR theory and its strongly coupled infrared modification which is likely 
to generate a stable ghost-free stage of cosmic acceleration and, perhaps, even simulate the DM effect 
on rotation curves. 

Of course, prospective nature of this model should not be exaggerated. Mainly, it challenges the 
ghost and fine-tuning problems in infrared modifications of Einstein theory, designed to simulate DE 
and DM mechanisms in the modern Universe. Without a proper extension this model most likely fails 
to describe early inflationary cosmology. In particular, it fails to generate the (quasi)de Sitter solution 
by the external cosmological constant A (vacuum energy of matter fields), because the relation p^ 
in this case becomes broken by a new scale X, a — —{A + 4i3)(l -I- A/A). Of course, this relation can 
be considered as the equation for A in terms of A and the stability analysis can be repeated, but the 
spirit of our approach avoiding the fine-tuning problem becomes lost. This property works in favor of 
a supersymmetric matter content of the Universe with a vanishing A. 

In contrast to a homogeneous vacuum energy source T^^ ex Xg^u, matter sources with a compact 
support present no difficulty and generate a stable perturbation theory on the (A)dS background, 
starting with the linearized potentials (j22p . This perturbation theory is applicable in the present 
day cosmology, and it should undergo tests on consistency of post-Newtonian corrections, on the 
magnitude of its DM phenomenon and the effect of the nonlocal trace term in (l22t . Moreover, the 
mechanism should be found, by which the model picks up a concrete scale |T^| ~ MpK of a crossover 
from the GR regime to cosmic acceleration - necessary element in realistic cosmologyo The last but 
not the least is the justification of the choice of ^ as an approximation for the effective action coming 
from some fundamental quantum gravity theory. These and other open issues that go beyond this 
paper will be discussed in [18]. 

In conclusion we mention that serendipity of ghost-free nonlocal gravity models Q satisfying the 
relations (|14p and ()18|) might not be exhausted by applications in cosmology. In particular, without the 
Riemann term in (|8|) (a = 0) they admit generic (not maximally symmetric) Einstein space solutions, 
R^J,u — ^g^J.v^ also with an arbitrary A, which might have implications for zero entropy black holes |21j 
and be an alternative to the conformal gravity model of 20 , whereas with a negative A they become 
a new testing ground for AdS/CFT correspondence perhaps promising other exciting consequences. 



^Disturbing property of perturbation theory in this phase could be the presence of poles ~ (l/n)ij in the vertices of 
(0 which could contribute too strong to graviton scattering. But these contributions vanish on shell /J^^ = provided, 
perhaps, that the Riemann term is forbidden in l[8}, a = 0. This question requires further studies. 

**Note that the existence of (A)dS solution with an arbitrary A is possible due to the fact that the purely gravitational 
action jTj transforms homogeneously under the global rescaling, g^i/ — > Ag^^, S[\gfj,i/\ 
mechanism can be based on the fact that matter fields break this property, cf. [?]• 
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